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ISSUE:  Several  methods  have  been  proposed 
for  determining  Manning’s  n  in  vegetated  chan¬ 
nels.  These  methods  recognize  that  die  physical 
charactedstics  of  the  vegetation  are  important 
factors  in  evaluating  flow  resistance.  However, 
minimal  research  has  been  conducted  to  quan¬ 
tify  the  density  and  drag  characteristics  of  vege¬ 
tation. 

RESEARCH  OBJECTIVES:  The  objectives 
of  the  research  presented  in  this  rqxirt  were  to 
develop  and  evaluate  methods  of  computing 
flow  resistance  in  vegetated  channels  and  flood- 
plains  and  to  characterize  the  impacts  of  vege¬ 
tation  for  flows  in  a  chaimel  typical  of  most 
temperate  floor^iains. 

SUMMARY:  A  field  study  was  performed  in  a 
vegetated  chaimel  located  near  Stillwater,  OK, 
as  a  part  of  the  development  of  a  comprehensive 
approach  to  predicting  resistance  to  flow  using 
the  Fischenich  equation.  The  charmel  was  char¬ 
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conditions:  (a)  without  leaves  on  trees  and 
shrubs  and  with  vegetative  debris  preset,  and 
(b)  with  leaves  on  trees  and  shrubs  and  without 
vegetative  debris. 
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used  to  develop  a  relation  between  the  coeffi¬ 
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drag  coefficient  associated  with  the  latter  condi¬ 
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specified  Re  was  10.2  times  greater  for  the  for¬ 
mer  condition  primarily  due  to  the  presence  of 
debris.  The  relations  for  dr^  can  be  input  into 
the  Fischenich  equation  to  solve  for  Maimmg’s 
n  for  channels  with  similar  vegetative  condi¬ 
tions.  Manning’s  n  was  predicted  within  an 
average  of  10  percent  for  four  streams  located 
in  Mississippi  and  Alabama  using  the  approach. 
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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  iinifs  of  measurement  used  in  this  report  can  be  converted  to  SI  umts  as 
follows: 


Multiply 

By 

To  Obtain 

acres 

4,046.873 

square  meters 

cubic  feet 

0.02831 

cubic  meters 

feet 

0.3048 

meters 

inches 

25.4 

millimeters 

miles  (U.S.  statute) 

1.609347 

kilometers 

pounds  (force  per 
square  foot) 

47.88026 

pascals 

pounds  (mass) 

0.4535924 

kilograms 

square  feet 

0.092903 

square  meters 

1  Introduction 


During  the  past  decade,  the  conventional  ‘‘flood  controF  ideology  has 
evolved  into  a  philosophy  of  “flood  management”  Streams  are  valuable  eco¬ 
nomic  assets  to  communities.  An  effective  flood  management  program  must 
consider  environmental,  recreational,  and  aesthetic  issues  in  addition  to  flood 
control.  Riparian  vegetation  is  an  integral  component  of  the  flood  channel.  Veg¬ 
etation  stabilizes  streambanks,  provides  shade  that  prevents  excessive  water  tem¬ 
perature  fluctuations,  supports  wildlife,  and  performs  an  essential  role  in  nutrient 
cycling  and  water  quality. 

Concurrent  with  the  benefits  provided  by  riparian  vegetation  are  the  issues 
pertaining  to  flood  hazard  mitigation.  Vegetation  increases  flow  resistance, 
which  has  an  inverse  effect  on  the  discharge  capacity  and  the  level  of  flood  pro¬ 
tection  provided  by  the  chaimel.  Several  methods  have  been  proposed  for  deter¬ 
mining  Manning’s  n  in  vegetated  channels.  These  mefliods  recognize  that  the 
physical  characteristics  and  growth  patterns  of  the  vegetation  are  important  fac¬ 
tors  in  evaluating  flow  resistance.  However,  little  guidance  is  available  to  quan¬ 
tify  the  physical  properties  of  the  vegetation  that  affect  flow  resistance.  A  coop¬ 
erative  study  between  the  U.S.  Army  Corps  of  Engineers  and  Colorado  State 
University  was  conducted  to  develop  a  method  of  characterizing  vegetation  for 
the  purpose  of  determining  resistance  to  flow  in  open  chaimels  or  floodplains. 

Background 

In  the  United  States  flow  resistance  is  customarily  expressed  in  terms  of  the 
resistance  coefficient  Manning’s  n.  The  resistance  coefficient  is  used  in 
Manning’s  Monomial  Equation  for  mean  flow  velocity  V  (Manning  1891) 

V  =  ^  R^JS^  (1) 

n  * 


where 

=  coefficient  equal  to  1.486  for  English  and  1  for  SI  units 
R  =  hydraulic  radius 
Sg  =  slope  of  energy  grade  line 
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The  primaiy  difficulty  in  using  Equation  1  is  <u;curately  estimating  the  value  of 
die  resistance  coefficient  Only  limited  guidance  is  available  for  selecting  resis¬ 
tance  coefficients,  particularly  in  vegetated  floodways. 


Fischenich  (1996)  presented  a  comprehensive  review  of  the  methods  available 
for  selecting  a  roughness  coefficient  for  vegetated  floodways  and  concluded  that 
the  existing  methods  do  not  permit  an  accurate  assessment  of  flow  depth  and 
velocity.  Based  on  conservation  of  momentum  principles,  he  suggested  that  an 
appropriate  relation  to  evaluate  flow  resistance  in  channels  with  nonsubmersed 
vegetation  is 


n 


KR 


2/3 


C,Veg, 


2g 


1/2 


(2) 


vdiere 

Q  =  coefficient  to  account  for  drag  characteristics  of  vegetation 
Vegj  =  vegetation  density 
g  =  gravitational  constant 

Vegetation  density  is  defined  as 


AL 


(3) 


where 

A,  =  area  of  vegetation  below  water  surface  projected  onto  a  plane 
perpendicular  to  direction  of  flow 

A  =  cross-sectional  area  of  flow 

L  —  characteristic  length,  such  as  a  unit  length  of  chamiel 

Equation  2  requires  knowledge  of  the  vegetation  density  and  drag  coefficient 
for  the  vegetation;  however,  litde  is  known  about  the  drag  characteristics  of  veg¬ 
etation.  The  vegetation  density  Vieg^can  be  directly  measured  using  a  horizontal 
point  frame  (Fhppin-Dudley  1997).  Figure  1  is  a  schematic  drawing  of  the  hori¬ 
zontal  point  frame.  The  frame  consists  of  two  vertical  poles  and  three  sets  (two 
each)  of  cross  arms,  1  m  long  and  inclined  at  an  angle  of  45  deg.  Each  cross  arm 
is  fabricated  with  equally  spaced  holes  that  are  aligned  with  the  holes  in  the  op¬ 
posite  cross  arm  through  which  a  pin  is  guided.  The  pin  is  0.3  m  long  and  sharp¬ 
ened  to  a  point  at  one  end.  Stakes  attached  to  the  lower  ends  of  the  poles  allow 
the  frame  to  be  freestanding  when  staked  into  the  ground. 

The  point  frame  is  first  staked  at  a  random  location  within  the  area  under 
investigation  (a  frame).  The  pin  is  slowly  guided  through  each  hole  (a  point)  in 
the  cross  arms  while  all  contacts  made  by  the  point  of  the  pin  with  vegetation  are 
recorded.  The  vegetation  density  Veg/v& 
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Number  of  Hits 
Number  of  Points 


J_ 

D 


(4) 


Veg.  = 


where  D  is  the  pin  length.  The  procedure  is  repeated  at  least  10  times,  or  until 
the  error  of  estimation  of  the  vegetation  density  for  the  delineated  area  as  deter¬ 
mined  by  the  running  mean  is  less  than  10  percent  Equation  4  applies  independ¬ 
ently  of  the  number  of  frames  or  points  per  frame  conducted.  For  example,  if 
200  hits  are  recorded  by  a  total  of  1,200  points  using  a  0.3-m  pin,  the  Veg^  is 
0.55/m  (200  hits/1,200  points/0.3  m). 

Study  Objective 

The  objective  of  the  study  was  to  develop  a  conq>rehensive  approach  using 
Equation  2  to  determine  the  Manning  roughness  coefficient  for  channels  with 
nonsubmersed  vegetation.  In  this  study,  the  coefficient  of  drag  was  calculated 
from  measured  values  of  Manning’s  n,  vegetation  density,  and  hydraulic  radius. 
Vegetation  and  flow  measurements  were  conducted  in  a  channel  located  at  the 
Agricultural  Research  Service  (ARS)  Outdoor  Hydraulic  LabOTatory  in 
Stillwater,  OK.  The  measurements  were  used  to  develop  a  relation  for  the  drag 
characteristics  of  the  vegetation  that  can  be  directly  input  into  Equation  2  for 

determining  Manning’s  n.  The  approach  was  then  substantiated  using  data  for 
four  streams  located  in  Mississippi  and  Alabama. 
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2  Field  Test  and  Application 


The  horizontal  point  frame  was  applied  to  measure  the  density  of  vegetation 
in  a  channel  located  at  the  ARS  Outdoor  Hydraulic  Laboratory  near  Stillwater, 
OK.  Vegetation  and  flow  measurennents  were  used  in  conjunction  with  Equa¬ 
tion  2  for  estimating  resistance  to  flow  (Manning’s  n)  for  channels  in  which  the 
vegetation  is  not  submersed.  In  Equation  2,  Manning’s  n  is  expressed  as  a  func¬ 
tion  of  the  flow  depth,  channel/flood^lain  geometry,  vegetation  density,  and  drag 
coefficient  In  practice,  flie  flow  depth  and  chaimel  geometry  are  usually  known. 
The  horizontal  point  frame  can  be  used  to  measure  the  vegetation  density  Vegj. 
Little  is  known  about  the  drag  characteristics  of  vegetation,  although  extensive 
information  is  available  concerning  the  coefficient  of  drag  for  common  objects 
and  shapes  (e.g.,  Blevins  1984). 

In  this  study,  the  horizontal  point  frame  was  used  in  conjunction  with  flow 
measurements  in  a  vegetated  channel  to  develop  a  relationship  for  the  drag  coef¬ 
ficient  of  vegetation  that  can  be  used  in  Equation  2  for  determining  Marming’s  n. 
Four  trips  were  taken  in  conjunction  with  the  testing — two  in  August  1995,  one 
each  in  November  1995  and  June  1996.  The  purpose  of  the  first  trip  in  August 
1995  was  reconnaissance  to  establish  locations  of  six  flow  measurement  cross 
sections.  During  the  August  1995  trip,  cross  sections  and  channel  slope  were 
surveyed,  plant  saiiqrles  were  collected,  and  preliminary  vegetation  measurements 
were  obtained.  Detailed  vegetation  measurements  were  conducted  in  June.  Flow 
measurements  were  conducted  in  both  November  1995  and  June  1996. 


Test  Channel 

History 

The  field  tests  were  conducted  at  the  ARS  Outdoor  Hydraulic  Laboratory 
located  approximately  10  km  west  of  Stillwater,  OK.  Figure  2  shows  the  test 
channel  in  November  1995  and  June  1996.  The  test  channel  was  constructed 
between  1942  and  1944  on  a  hillside  composed  of  heavy  clay  soil  at  the  base  of 
the  Carl  Blackwell  Dam.  Water  is  conveyed  from  the  lake  through  five  siphons, 
each  having  a  maximum  capacity  ranging  between  0.57  and  0.85  m^/s,  depending 
on  the  water  level  of  the  reservoir.  The  siphon  flow  rate  is  measured  using  an 
ogee  crest  weir  structure  at  the  outlet.  Downstream  from  the  weir,  the  water 
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a.  November  1995 


b.  June  1 996 

Rgure  2.  Test  channel  in  Stillwater,  OK 


flows  by  gravity  through  a  main  distribution  channel  and  enters  the  test  channel 
through  a  3-m  sluice  gate.  At  the  terminus  of  the  test  channel,  the  water  is  re¬ 
leased  into  another  chaimel  that  discharges  into  Stillwater  Creek. 
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The  original  channel  was  117  m  long  and  divided  into  three  reaches  by  slope. 
The  upper  reach  was  originally  graded  at  a  3-percent  slope  for  41m.  Hie 
47-m-long  middle  reach  was  constructed  at  a  6-percent  grade  and  the  30-m-long 
lower  reach  was  graded  at  a  4-percent  slope.  The  cross  section  was  V-shaped, 
0.6  m  deep,  and  12  m  wide  wife  1:10  side  slopes  (Gwiim  and  Ree  1980). 

During  the  period  from  1944  through  1950,  the  chaimel  was  maintained  by 
mowing  and  seeding  with  various  types  of  grasses  in  association  with  a  research 
program  designed  to  determine  the  stability  and  capacity  of  a  grass-lined  chan¬ 
nel.  The  results  of  this  research  were  the  ARS  n-VR  method  for  determining 
Manning’s  n  for  grass-lined  chaimels  (U.S.  Department  of  Agriculture  (USD A) 
1954).  The  chaimel  has  not  been  maintained  since  1950,  and  volunteer  vegeta¬ 
tion  (i.e.,  trees,  shrubs,  and  grasses,  etc.)  has  been  permitted  to  grow  without 
interference. 

Plant  composition 

Since  1951,  the  chaimel  has  become  overgrown  with  vegetation  typical  of 
dense  woods  and  thickets  found  in  northwestern  Oklahoma.  Nineteen  plant 
specimens  were  collected  during  August  1995.  They  consisted  of  eight  tree,  six 
grass  and  forb,  one  shrub,  and  two  woody  vine  species  (Tables  1  and  2). 

The  tree  species,  percent  composition,  and  average  diameter  of  each  tree  spe¬ 
cies  are  presented  in  Table  1.  The  average  diameter  of  all  trees  measured  was 
10  cm.  Only  trees  having  diameters  greater  than  2J5  cm  were  measured.  Dog¬ 
wood  (Comus  drummondii),  elm  (Ulmus  nubra),  and  hackbeny  {Celtis 


Table  1 

Composition  and  Average  Diameter  of  Trees  in  Stillwater,  OK,  Test 
Channel  


Common  Name 

Scientific  Name 

Composition,  % 

Average  Diameter,  m 

Dogwood 

Comus  drummondii 

33 

0.05 

Elm’ 

Ulmus  nubra 

21 

0.19 

1  Hackbeny 

Celtis  tenurfotia 

Red  Cedar 

Juniperus  virginiana 

15 

0.07 

Honey  Looet 

Gleditsia  triacanthos 

11 

0.17  1 

Unkno¥m  (dead) 

Unknown 

8 

0.04  1 

Pecan 

Carya  pecan 

5 

0.11  1 

Upland  Sumac^ 

Rhus  copal fina 

5 

0.04 

Smooth  Upland  Sumac 

Rhus  glabra 

Willow 

Salixspp. 

2 

0.19 

’  No  distinction  made  between  elm  and  hackbeny. 

^  No  distinction  made  between  upland  and  smooth  upland  sumac. 
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b  Table  2 


Composition  of  Grasses  and  Foi1>s,  Shrubs,  and  Woody  Vines  in 
Stiilwater,  OK,  Test  Channei 


Plant  Group 

Common  Name 

Scientific  Name 

Composition,  % 

Grasses  and  Forbs 

Verrain 

Prairie  Sage 

Poverty  Oak  Grass 
Sedge 

Canada  Goldenrod 

Rice  Cut-Grass 

Verbena  scabra 

Artemisia  eudoviciana 
Danthonia  spicata 

Carexspp, 

Sotidago  canadersis 

Leersia  virginica 

50  1 

Shrubs 

Huckfeberry 

Vaccinium  pallidum 

36  1 

Woody  Vines 

Carofina  Moonseed 
Poison  Ivy 

Cocculus  carofmus 
Toxicodendron  rydbergu 

14 

temdfolia)  comprised  more  than  50  percent  of  the  trees  in  the  test  channel.  These 
species  were  present  throu^out  the  entire  length  of  the  channel.  Red  cedar 
(Jumperus  virginiana),  comprising  15  percent  of  the  trees,  was  most  prevalent  in 
the  lower  reach  of  the  charmel.  The  red  cedar  exposed  to  the  sun  on  the  charmel 
margins  speared  healthy.  Cedar  trees  inside  the  canopy  were  devoid  of  leaves 
and  often  dead,  ^parently  due  to  a  lack  of  sunlight  or  flooding.  Sumac  (Rhus 
glabra  and  Rhus  copaUind)  and  honey  locust  (Gleditsia  triacanthos)  comprised 
5  and  1 1  percent  of  the  trees,  respectively,  and  were  common  in  the  upper  and 
middle  reaches  of  the  charmel.  Pecan  (Carya  iUinoensis)  and  wiUow  (Salix  spp.) 
comprised  5  and  2  percent  of  die  trees  in  the  charmel,  respectively,  and  were  only 
encountered  in  the  lower  reach  of  the  channel. 

The  term  composition  has  a  different  meaning  in  Table  1  than  in  Table  2  be¬ 
cause  the  tree  group  was  evaluated  separately  from  fee  other  three  plant  groups. 

In  Table  2,  fee  composition  is  fee  percentage  of  fee  stand  comprised  by  each 
plant  group.  That  is,  50  percent  of  fee  vegetation  encountered  (excluding  trees) 
was  grasses  and  forbs,  36  percent  was  shrub,  and  14  percent  was  woody  vine. 
Grasses  and  forbs  were  found  feroughout  fee  diannel,  particularly  in  areas  ex¬ 
posed  to  sunlight  such  as  fee  charmeTs  edge,  but  were  concentrated  in  fee  area  of 
fee  grade  break  between  the  upper  and  middle  reacl^  and  near  fee  charmeTs 
end.  None  of  fee  original  grasses  planted  in  fee  charmel  as  a  part  of  the  research 
conducted  in  fee  1940s  were  observed.  Ihe  only  shrub  encountered  was  fee 
huckleberry  (Vaccinium  pallidum),  which  was  prevalent  feroughout  fee  channftl 
The  two  woody  vines,  poison  ivy  (Toxicodendron  rydbergii)  and  Carolina  mcxin- 
seed  (Cocculus  carolinus),  were  abundant  in  fee  lower  reach. 

Channel  survey 

Six  cross  sections  were  established  and  surveyed  in  August  1995  (Figure  3). 
The  surveyed  charmel  profile  and  a  typical  cross  section  of  fee  test  charmel  are 
shown  on  Figures  4  and  5,  respectively.  The  charmel  profile  and  geometry  are 
similar  to  fee  original  shape.  The  1 17-m-long  charmel  can  be  divided  into  three 
reaches  according  to  slope.  Cross  sections  1  through  4  are  ^proximately  12  m 
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Rgure  4.  Profile  of  channel  showing  lower,  middle,  and  upper  reaches  and  six 
cross  sections 


wide,  0.6  m  deep,  with  1:10  side  slopes.  The  channel  gradually  widens  to  ap- 
proximately  18m  near  sections  5  and  6. 

A  strip  gully,  approximately  0.3  to  0.6  m  wide,  forms  a  well-defined  thalweg 
that  roughly  follows  the  center  of  the  chaimel  (Figure  3).  The  gully  meanders 
around  trees  and  large  roots  and  was  likely  formed  by  runoff  during  rainstorms  or 
leakage  under  the  sluice  gate  from  the  main  distribution  chaimel. 
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Rgure  5.  Typical  cross  section,  Stillwater.  OK,  test  channel 

Vegetation  Density  Measurements 


Vegetation  density  measurements  were  conducted  in  the  test  channel  in  Au¬ 
gust  1995  and  June  1996.  During  die  initial  assessment,  measurements  were 
taken  at  each  of  the  six  cross  sections,  and  tree  diameters  were  recorded  at  ran¬ 
dom  locations.  These  data  were  used  to  assess  die  general  character  of  vegeta¬ 
tion  in  die  channel.  A  more  intense  examination  of  the  vegetation  was  conducted 
during  June  1996;  these  latter  measurements  served  as  the  basis  for  the  vegeta¬ 
tion  analysis  presented. 

Vegetation  measurements  were  taken  up  to  a  height  of  0.6  m  above  the 
thalweg  because  the  maximum  depth  of  flow  in  the  test  channel  was  not  expected 
to  exceed  this  height  Measurements  were  conducted  using  two  techniques: 

(a)  the  horizontal  point  frame  and  (b)  direct  measurement  of  tree  trunks.  Grass 
and  forb,  shrub,  and  vine  plant  groups  were  measured  using  the  horizontal  point 
frame.  Trees  were  measiured  by  both  the  horizontal  point  frame  method  and  di¬ 
rect  measurement  of  diameter.  All  leaves,  stems,  and  tree  parts  and  trunks  less 
than  2.5  cm  in  diameter  were  included  in  the  point  frame  measurements.  Tree 
trunks  greater  than  2.5  cm  in  diameter  were  measured  directly  in  November  1995 
(when  leaves  on  the  trees  and  shrubs  were  absent)  to  determine  the  Vegj. 
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Horizontal  point  frame  measurements 

Horizontal  point  frame  measurements  were  taken  using  the  line  transect 
method.  The  line  transect  method  is  a  systematic  sampling  technique  that  con¬ 
sists  of  establishing  a  line  across  the  channel  and  placing  the  point  frame  at  spe¬ 
cific  intervals  along  the  line.  Thirteen  line  transects  were  established  (Rgure  6). 
The  line  transects  were  marked  by  a  cloth  surveying  tape  extended  across  the 
channel  at  a  height  of  0.6  m  above  the  thalweg.  The  length  of  the  transects  var¬ 
ied  from  9  to  14  m  according  to  the  channel  cross-sectional  characteristics  at 
each  transect  (Table  3). 


20  40  60  80  100  120 

LoBgItudilud  Distance  (m) 

Top  of  Bank  X  Frame  Mcasnrement  limit  of  Line  Transect 

Line  Transect  - Cross -section 


Figure  6.  Plan  wew  of  channel  depicting  six  cross  sections  and  1 3  line 
transects  along  which  vegetation  density  measurements  were 
obtained  using  horizontal  point  framne 


Four  to  six  point  frame  measurements  2  m  apart  were  taken  along  each  tran¬ 
sect,  and  Table  3  summarizes  the  data  for  the  13-line  transects.  Measurements 
were  talfp.n  with  the  lowest  cross  arm  of  the  horizontal  point  frame  (Figure  1) 
since  die  maYimnin  depth  of  flow  was  expected  to  be  0.6  m.  The  frame  was 
placed  adjacent  to  the  line,  and  points  were  inserted  up  to  die  height  of  the  t^. 
The  number  of  points  used  per  frame  depended  on  the  height  of  the  tape  above 
the  charmel  bottom  at  each  location,  and  varied  from  a  minimum  of  4  points/ 
frame  near  the  edges  of  the  charmel  to  a  maximum  of  20  points/frame  at  the 
thalweg. 

Out  of  a  total  of  906  points  (67  frames),  158  hits  were  recorded.  Sixty-eight 
hits  were  made  with  grass  or  forb,  19  with  vines,  48  with  shrub,  and  23  with  trees 
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Table  3 

Summary  of  Horizontal  Point  Frame  Measurements  Conducted  In  Stillwater.  OK.  Test  Channel 


(excluding  contacts  with  tree  trunks  greater  than  2.5  cm  in  diameter).  ITie  den¬ 
sity  of  the  vegetation  measured  with  the  point  fiame  as  calculated  from  Equa¬ 
tion  4  is  0.581/m  (158  hits/906  points/0.3  m). 

Tree  measurements 

Belt  transects  were  established  at  line  transects  2  through  13  for  the  purpose 
of  determining  the  contribution  of  tree  trunks  greater  than  2.5  cm  in  diameter.  A 
belt  transect  is  a  rectangular-shaped  or  strip  quadrat  The  belt  transects  in  this 
study  were  oriented  with  die  long  dimension  perpendicular  to  the  flow  direction. 
The  dimensions  were  9  by  3  m  (1.5  m  upstream  and  downstream  of  the  line).  A 
belt  transect  was  not  established  in  conjunction  with  the  first  line  transect  due  to 
its  close  proximity  to  the  channel  entrance  gate.  The  combined  width  of  the 
12  belt  transects  (36  m)  represented  over  30  percent  of  the  total  charmel  length 
(117  m).  The  tree  species,  location  with  the  belt  transect,  basal  diameter,  and 
diameter  at  0.6-m  height  were  recorded.  A  total  of  126  trees  having  a  trunk  di¬ 
ameter  greater  than  2.5  cm  were  identified  within  the  12  transects.  Hgure  7  pres¬ 
ents  a  plan  view  of  the  test  charmel  showing  the  location  of  the  12  belt  transects 
(labeled  2  through  13)  and  trees  in  which  the  diameters  were  directly  measured. 


Rgure  7.  Plan  view  of  test  channel  depicting  six  cross  sections  and  12  belt 
transects  in  w4iich  tree  trunks  >2.5  rxn  were  directly  measured 

Table  4  contains  a  surmnary  of  the  calculations  to  determine  the  contribution 
of  the  trees  to  die  total  Vegj .  The  Vegj  of  the  tree  trunks  with  greater  than 
2.5  cm  in  diameter  was  calculated  as 
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Table  4 

Summary  of  Veg^  Calculations  for  Tree  Trunks  Greater  Than 
2.5  cm  in  Diameter 


Belt  Transect  No. 

lA 

Ag 

Veg^ 

1/m 

2 

0.611 

4.279 

0.048 

3 

0.634 

3,337 

0.063 

4 

0.930 

3.770 

0.082 

5 

0.395 

3.682 

0.036 

6 

0.258 

4.020 

0.021 

7 

0.055 

3.214 

0.006 

8 

0344 

3.590 

0.051 

10 

0.652 

2.665 

0.082 

11 

0.855 

3.907 

0.073 

12 

0.232 

4.553 

0.017 

13 

0.055 

3.958 

0.005 

Vegj 


A^L 


(5) 


where 

Aj  =  projected  area  of  an  individual  tree  below  water  surface 
Ag  =  area  of  flow  represented  by  belt  transect 
L  =  length  of  belt  transect,  3  m 

The  value  of  A,  was  calculated  by 


=  yi^i 


(6) 


where 

=  depth  of  flow  at  individual  tree 
=  mean  of  tree  basal  diameter  and  diameter  at  0.6  m 

The  flow  depth  at  each  tree  y,-  depends  on  the  horizontal  station  of  the  tree  within 
the  transect  because  the  channel  cross  section  is  generally  triangular  in  shape. 
The  value  of  y,  was  determined  from  the  vertical  distance  between  the  tape  that 
marked  the  line  transect  (strung  at  a  height  of  0.6  m  above  the  thalweg)  and  the 
base  of  the  tree.  The  characteristic  flow  area  Ag  was  calculated  based  on  the 
channel  geometry  at  the  line  transect  and  the  width  of  the  belt  transect  (9  m). 
Table  4  lists  the  total  area  of  all  trees,  flow  area,  and  vegetation  density 
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associated  with  each  belt  transect  The  mean  Veg^  of  the  entire  test  channel  due 
to  tree  trunks  greater  than  2.5  cm  is  0.046/m.  Summing  the  tree  and  horizontal 
point  measurements  yields  a  mean  for  the  entire  channel  of 0.627/m 

(0.046/m  +  0.581/m). 

The  values  of  Veg^  for  November  were  determined  based  on  measurements 
conducted  in  June.  Tte  tree  trunk  measiuements  were  used  directly  because  the 
contribution  of  the  tree  tnmks  to  the  value  of  Vcg^does  not  change  between 
seasons.  The  primary  difference  between  the  condition  of  the  vegetation  in  No¬ 
vember  and  June  was  that  the  tree  and  shrub  leaves  were  not  present.  Therefore, 
all  contacts  with  leaves,  including  grass,  were  excluded  from  the  horizontal  point 
frame  measurements.  Grass  was  present  in  the  channel  during  the  November 
test;  however,  it  was  not  included  in  the  Vegj  parameter  because  the  density  of 
the  grass  was  not  known.  The  effect  of  excluding  the  grass  on  the  analysis  of 
resistance  are  discussed  in  more  detail  in  the  section  on  assmnptions  of  the 
methodology.  The  Vcg^  used  in  the  analysis  for  November  was  0.123/m  based 
on  the  tree  trunk  measurements  and  contacts  made  witfi  stems  using -the  point 
frame  (0.077/m  stems  +  0.046/m  tree  trunks). 

Assessment  of  vegetation  density 

FUppin-Dudley  et  al.  (in  preparation)  recommended  three  st^s  for  field  ^)pli- 
cation  of  die  horizontal  point  frame:  (a)  delineation  of  area,  (b)  location  of 
measurements,  and  (c)  decision  of  san^le  size.  The  steps  were  applied  to  the 
vegetation  measurements  conducted  in  the  test  channel  as  follows: 

CL  Delineation  of  Area  -  The  area  in  which  the  measurements  are  conducted 
should  correspond  to  the  length  of  channel  in  which  the  roughness  charac¬ 
teristics  are  to  be  evaluated.  In  this  study,  the  measurements  were  con¬ 
ducted  throughout  the  entire  length  of  the  test  channel.  An  average  vege¬ 
tation  density  was  used  because  die  vegetation  (resistance)  characteristics 
were  similar  throughout  the  channel. 

b.  Location  of  Measurements  -  Vegetation  measurements  were  conducted 
based  on  systematic  sampling  methods  using  line  transects  for  hc^izontal 
point  frame  data  and  belt  transects  for  tree  data.  Measurements  should  be 
obtained  by  random  sampling  in  the  area  under  investigation  to  minimize 
the  effects  of  trends  in  the  vegetation;  therefore,  the  data  were  randomized 
prior  to  analysis  of  sarr^le  size. 

c.  Sample  Size  -  A  plot  of  the  running  mean  was  used  to  assess  die  accuracy 
of  die  mean  vegetation  density.  Uie  method  consists  of  plotting  succes¬ 
sive  or  “rurming”  means  against  the  number  of  measurements.  Sampling 
should  be  continued  until  the  percentage  difference  between  the  previous 
mean  and  the  new  mean  is  5  percent  or  less.  The  minimum  number  of 
frames  should  correspond  to  that  required  to  obtain  either  100  hits  with 
vegetation  or  10  frames. 
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Plots  of  the  running  mean  were  prepared  for  the  randomized  horizontal 
point  fiame  and  tree  data  (Figure  8).  The  analysis  of  the  data  indicated 
that  the  difference  in  cumulative  mean  values  was  2  percent  after 
30  frames.  However,  the  number  of  hits  with  vegetation  for  30  frames  was 
56,  which  is  less  than  the  specified  mininwiTn  of  100.  The  number  of 
fiames  required  to  obtain  100  hits  with  vegetation  was  42.  In  this  study, 
67  fiames  were  conducted,  which  is  greater  than  the  Tninimnm  of  42. 

A  mean  vegetation  density  was  plotted  for  the  tree  data  based  on  the  mea¬ 
surements  conducted  for  each  transect  The  difference  in  cumulative  mean 
values  was  4  percent  after  eight  belt  transects.  The  number  of  belt  tran¬ 
sects  required  to  obtain  a  count  of  100  trees  was  1 1.  The  actual  number  of 
belt  transects  analyzed  was  12. 

Flow  Depth  and  Velocity  Measurements 

How  depth  and  velocity  measurements  were  conducted  in  November  1995 
and  June  1996.  It  was  desired  to  determine  the  flow  resistance  characteristics  of 
the  channel  in  its  existing  condition  in  November  and  then  again  after  the  re¬ 
moval  of  vegetative  debris  and  growth  of  green  foliage  in  June.  In  both 
November  and  June,  flow  depth  and  velocity  measurements  were  conducted  for 
nine  test  discharges  at  each  of  the  six  cross  sections.  Each  test  discharge  was 
preliminarily  established  based  on  static  head  measurements  for  a  calibrated  ogee 
crest  weir  located  upstream  of  the  channel  entrance.  The  ninth  discharge  was 
less  during  June  testing  due  to  lower  water  levels  in  the  water  supply  lake. 

The  flow  depth  and  velocity  measurements  were  conducted  according  to  the 
procedures  presented  in  the  National  Handbook  of  Recommended  Methods  for 
Water-Data  Acquisition^  U.S.  Geological  Survey  (USGS)  (1977).  First,  a  tag  line 
(survey  tape)  was  extended  across  the  water  surface  to  marie  the  cross  section. 

The  flow  depth  was  measured  using  a  wading  rod  placed  at  even  increments 
along  the  tag  line.  Tlie  increments  ranged  from  0.3  to  0.6  m  depending  on  the 
width  of  the  water  surface.  At  each  location,  the  velocity  was  measured  using  a 
Marsh-McBimey  velocity  meter.  Velocity  measurements  were  conducted  using 
the  six-tenths  method  when  the  flow  depth  was  less  than  0.3  m  or  two-point 
method  when  the  flow  depth  equaled  or  exceeded  0.3  m.  In  the  six-tenths 
method,  the  velocity  is  recorded  at  a  point  located  six-tenths  of  the  total  depth 
below  the  water  surface.  In  the  two-point  method,  the  velocity  is  recorded  at 
two-tenths  and  eight-tenths  of  the  total  depth  below  the  water  surface.  The  aver¬ 
age  of  die  two-point  measurements  is  used  to  compute  discharge.  The  average 
water  tenqierature  was  13  in  November  and  25  in  Jime. 

Calculation  of  Discharge 
and  Manning’s  n 


The  discharge  was  calculated  using  the  velocity-area  method  and  the  midsec- 
tion  method  for  computing  cross-sectional  area  (USGS  1977).  In  the  midsection 
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Mean  Veg^  (1/m)  Mean  Veg^  (1/m) 


Number  of  Frames 

(a)  Hoii2ontal  point  frame  measurements 


(b)  Tree  measurements 


Rgure  8.  Plots  of  mean  mnning  area 
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method,  the  meastired  velocity  and  flow  depth  are  assumed  to  represent  the  mean 
values  for  a  rectangular  area  extending  laterally  one-half  the  distance  to  the  adja¬ 
cent  measiu'ement  stations.  The  total  discharge  is  the  summation  of  the  products 
of  the  partial  areas  (a^  and  velocities  (v^)  for  each  cross  section. 


Manning’s  n  was  calculated  based  on  the  slope  of  the  energy  grade  line,  be¬ 
tween  two  adjacent  cross  sections 


Lwse+hk. 


(7) 


where 

Awse  =  difference  in  water  surface  elevations 
Ah^  =  difference  in  velocity  head 

L  =  length  of  reach 


The  velocity  head  at  each  section  was  calculated  as 


h 


V 


2g 


(8) 


where  a  is  the  energy  coefficient  The  value  of  a  was  computed  as 

^  (Sv/fl.)(Sv/ 


(9) 


From  Equation  7  and  Manning’s  equation  (Equation  1),  Manning’s  n  for  a  chan¬ 
nel  reach  containing  two  or  more  cross  sections  is  expressed  as 


n 


J_ 

Q 


(wsej-wse^)+(h^-h^) 


^2, 


L: 


V-I)j 


where  Z  is  a  section  factor 


2/3 


(10) 


(11) 


and  j  is  the  number  of  cross  sections  (which  are  designated  as  1, 2, 3, ... 
with  tile  jlii  cross  section  being  the  farthest  upstream).  Equation  10  is  used  to 
calculate  Manning’s  n  for  a  specified  channel  reach  (e.g.,  between  cross  sections 
or  for  multiple  reaches  between  cross  sections).  Manning’s  n  was  calculated 
for  the  entire  channel  and  five  reaches  shown  on  Figure  3.  Tables  5  and  6 
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Table  5 

Average  Hydraulic  Properties,  Stiilwatei 

r,  OK,  Te 

!St  Chan 

nei,  Nov 

ember  1 

995  .] 

Row 

BM 

Reach 

L 

HH 

R 

H 

No. 

No. 

m 

m 

a 

n  1 

1 

0.10 

1 

6 

0.16 

0.36 

0.08 

0.28 

2.55 

0.039 

0.128  1 

7 

0.20 

0.36 

0.09 

0.27 

3.22 

0,031 

0.123  1 

19 

0.21 

0.35 

0.11 

0.28 

4.43 

0.040 

0.148  I 

52 

0.21 

0.42 

0.10 

0.27 

3.86 

0.055 

0.201  1 

10 

0.24 

0.49 

0.12 

0.23 

2.70 

0.035 

0.209 

2 

0.33 

1 

6 

0.25 

0.86 

n 

0.41 

3.73 

0.038 

0.145 

7 

0.32 

0.91 

BH 

0.36 

3.33 

0.024 

0.126 

19 

0.35 

0.93 

0.33 

3.41 

0.041 

0.181 

52 

0.30 

1.03 

0.17 

0.33 

3.07 

0.056 

0.225 

10 

0.35 

1.20 

0.17 

0.27 

4.99 

0.029 

0.193 

3 

0.64 

1 

6 

0.32 

1.26 

0.19 

0S4 

2.93 

0.047 

0.138  1 

7 

0.40 

1.52 

0.21 

0.48 

3.40 

0.016 

0.102 

19 

0.47 

1.58 

0.23 

0.45 

3.07 

0.041 

0.183 

52 

0.40 

1.24 

0.19 

0.53 

2.87 

0.057 

0.153 

10 

0.40 

1.13 

0.18 

0.53 

3.15 

0.032 

0.099 

4 

0.85 

1 

0.37 

1.46 

0.20 

0.57 

3.15 

0.049 

0.130 

0.43 

1.77 

0.24 

0.51 

2.85 

0.021 

0.114 

0.49 

1.77 

0.26 

0.49 

3.03 

0.039 

0-162 

0.43 

1.45 

0.21 

0.55 

3,02 

0.057 

0.145 

0.44 

1.77 

0.21 

0.50 

3.58 

0.030 

0-127 

5 

1.34 

1 

6 

0.43 

1.90 

0.24 

0.72 

2.57 

0.040 

0.107 

7 

0.52 

2.34 

0.27 

0.61 

2,50 

0.022 

0.105 

19 

035 

2.28 

0.26 

0.62 

2.49 

0.042 

0.137 

52 

0.49 

1.92 

0.23 

0.72 

2.79 

0.056 

0.125 

10 

0.53 

2.39 

0.26 

0.57 

4.02 

0.033 

0.129 

6 

1.78 

1 

6 

0.49 

2.26 

0.80 

1.86 

0.047 

0.111 

7 

0.52 

2.70 

0.67 

2.27 

0.030 

0.109  1 

19 

0.53 

2.65 

0.26 

0.70 

2,10 

0.041 

0.119  1 

52 

032 

2.36 

0.24 

0.78 

2-42 

0.054 

0.119  1 

1^1 

10 

0.55 

2.54 

0.26 

0.70 

2.74 

0.042 

0.119  1 

7 

2.28 

6 

030 

2.68 

0.28 

0,87 

1.68 

0.050 

0.112 

7 

0.55 

2.87 

0.29 

0.81 

1.86 

0.022 

0.081 

19 

0.56 

2.82 

0.28 

0.81 

2.02 

0.043 

0.109 

52 

0.56 

2.87 

0.30 

0.82 

2.67 

0.053 

0.127 

10 

0.64 

3.10 

0.32 

0.75 

2.90 

0.040 

0.127 

8 

3.07 

1 

6 

0.56 

3.10 

0.30 

1.01 

1.88 

0.034 

0.084  1 

7 

0.61 

3.33 

0.32 

0.94 

1.76 

0.039 

0.098  1 

19 

0.62 

3.22 

0.30 

0.97 

1.^ 

0.037 

0.089 

52 

0.62 

3.15 

0.30 

0.98 

2.29 

0.055 

0.106 

10 

0.66 

3.48 

0.32 

0.88 

2.73 

0.036 

0.101 

9 

3.97 

1 

6 

0.61 

3.63 

0.33 

1.09 

1.73 

0.060 

0.107 

7 

0.64 

4.00 

0.34 

1,00 

1-55 

0.023 

0.074 

19 

0.67 

3.80 

0.33 

1.10 

1.54 

0.040 

0.090 

52 

0.64 

3.73 

0.33 

1.11 

I  2.30 

0.055 

0.102  1 

5  J 

10 

0.69 

4.15 

0.34 

0.96 

2.79 

0.036 

0.096  H 
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Table  I 
Averac 

5 

Hydra 

ulic  Pro| 

perties,  1 

Stillwate 

r,  OK,  Test  Channel,  June  1 996 

Flow 

Q 

Reach 

L 

BH 

lilH 

B 

PU 

III^H 

\mm 

No. 

mVs 

No. 

m 

HH 

m 

1 

0.11 

1 

6 

0.18 

0.41 

0.09 

0.26 

2.88 

0.045 

— 

0.150 

7 

0.18 

0.32 

0.08 

0.35 

2.70 

0.034 

0.099 

19 

0.23 

0.29 

0.09 

0.42 

2.95 

0.035 

0.096 

52 

0.23 

0.30 

0.09 

0.39 

3.15 

0.057 

0.125 

\wam 

0.20 

0.29 

0.09 

0.38 

2.65 

0.042 

0.108 

2 

0.34 

6 

0.26 

0.74 

0.13 

0.42 

2.23 

0.035 

0.105 

7 

0.27 

0.73 

0.13 

0.42 

2.38 

0.045 

0.117 

19 

0.32 

0.71 

0.14 

0.50 

2.51 

0.035 

0.106 

52 

0.30 

0.62 

0.14 

0.58 

2.14 

0.056 

0.117 

1 

wsmm 

0.27 

0.55 

0.13 

0.66 

2.07 

0.034 

0.175 

fs 

0.66 

6 

0.34 

1.16 

0.18 

0.58 

1.82 

0.054 

0.130 

n 

7 

0.34 

1.18 

0.19 

0.55 

2.12 

0.029 

0.100 

fl 

19 

0.35 

1.12 

0.19 

0.66 

2.23 

0.039 

0.108 

1 

52 

0.35 

0.89 

0.17 

0.82 

1.81 

0.055 

0.098 

■■I 

0.37 

0.87 

0.16 

0.72 

2.06 

0.044 

0.081 

0.93 

1 

6 

0.38 

1.41 

0.20 

0.62 

2.11 

0.052 

0.118 

1 

7 

0.41 

1.51 

0.22 

0.59 

1.72 

0.026 

0.094 

19 

0.41 

1.42 

0.21 

0.72 

1.92 

0.041 

0.105 

52 

0.38 

1.11 

0.18 

0.90 

2.05 

0.055 

0.087 

10 

0.41 

1.13 

0.18 

0.86 

2.31 

0.039 

0.076 

5 

1.39 

1 

6 

0.44 

1.79 

0.23 

0.80 

1.72 

0.051 

0.109 

7 

0.43 

1.85 

0.25 

0.76 

1.48 

0.039 

0.104 

19 

0.43 

1.75 

0.24 

0.90 

1.63 

0.036 

0.088 

52 

0.40 

1.31 

0.19 

1.05 

1.72 

0.057 

0.075 

10 

0.41 

1.39 

0.20 

0.94 

1.93 

0.034 

0.062 

6 

1.86 

1 

6 

0.47 

2.11 

0.87 

1  1.93 

0.051 

0.101 

7 

0.46 

2.08 

0.89 

1.42 

0.039 

0.092 

19 

0.46 

1.95 

0.94 

1.68 

0.036 

0.077 

52 

0.43 

1.61 

0.22 

1.18 

1.75 

0.056 

0.075 

BHH 

10 

0.44 

1.77 

0.23 

1.14 

1.82 

0.037 

0.069 

7 

228 

H 

6 

0.49 

2.48 

0.27 

0.94 

1.86 

0.056 

0.107 

7 

0.47 

2.41 

0.28 

0.95 

1.41 

0.036 

0.085 

19 

0.47 

2.26 

0.26 

1.04 

1.42 

0.038 

0.077 

52 

0.44 

1.92 

022 

1.25 

1.59 

0.055 

0.072 

10 

0.52 

2.07 

028 

1.12 

1.78 

0.032 

0.061 

8 

2.97 

1 

6 

0.53 

3.11 

0.30 

1.02 

1.79 

0.049 

0.102 

7 

0.56 

3.02 

0.30 

1.00 

1.55 

0.029 

0.077 

19 

0.62 

2.87 

0.28 

1.04 

1.56 

0.035 

0.075 

52 

0.56 

2.49 

0.25 

1.18 

1.80 

0.057 

0.080 

^mi 

10 

0.56 

2.66 

0.28 

1.10 

1-92 

0.033 

0.070 

9 

3.42 

1 

6 

0.59 

3.47 

0.32 

0.99 

1.80 

0.052 

0.108 

7 

0.58 

3.43 

0.32 

0.98 

1.58 

0.038 

0.092 

19 

0.58 

3.17 

0.30 

1.10 

1.58 

0.037 

0.078 

52 

0.56 

2.83 

0.27 

1.21 

1.95 

0.055 

0.081 

■■■11 

10 

0.56 

3.04 

0.30 

1.13 

2.02 

0.046 

0.085 
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summarize  the  values  of  Manning’s  n  and  other  average  hydraulic  properties  for 
the  five  reaches.  The  Manning’s  n  values  for  the  entire  channel  based  on  the 
six  cross  sections  are  listed  in  Table  7. 


Table? 

Discharge  and  Manning’s  n  Values  for  Stillwater,  OK,  Test 

Channel 

Flow  No. 

November^ 

June^ 

Difference,  % 

O  E 

n 

Q 

mVs 

n 

1 

0.10 

0.176 

0.11 

0.117 

51 

2 

0.33 

0-199 

0.34 

0,109 

83 

3 

0.64 

0.147 

0.66 

0.099 

48 

4 

0.85 

0.143 

0.93 

0.090 

59  1 

5 

1.3 

0-125 

1.4 

0.078 

60  1 

6 

1.8 

0.118 

1.9 

0.077 

53 

7 

2.3 

0.119 

2.3 

0-074 

61 

8 

3.1 

0.100 

3.0 

0.079 

27 

9 

4.0 

0.098 

3.4 

0.083 

18 

'  Debris  preserrt  and  leaves  absent.  | 

^  Debris  absent  and  leaves  present.  1 

Discussion  of  Results 

Evaluation  of  Manning’s  n 

Table  8  presents  the  percentage  difference  in  Manning’s  n  values  for  the 
five  reaches  and  the  nine  flow  tests  conducted  in  both  November  and  June.  The 
values  are  expressed  as  a  percentage  of  the  June  value.  The  two  primary  differ¬ 
ences  in  the  channel  conditions  between  the  two  periods  are  that,  in  November, 
the  trees  and  shrubs  were  devoid  of  leaves,  and  vegetative  debris,  such  as  dry 
leaves  and  wood,  were  present  The  values  in  Table  8  reflect  both  the  influence 
of  the  debris,  which  would  increase  the  value  of  Maiming’s  n,  and  the  lack  of 
leaves,  which  would  decrease  Manning’s  n  relative  to  the  June  value. 

It  is  apparent  from  Tables  7  and  8  that  the  presence  of  debris  had  significant 
influence  on  the  flow  resistance  and  measured  value  of  Manning’s  n.  Table  7 
indicates  an  increase  in  the  average  percent  difference  in  Manning’s  n  from  up¬ 
stream  to  downstream  (reach  1  to  5).  The  average  percent  difference  tends  to 
decrease  with  increasing  discharge  particularly  for  flow  numbers  7  through  9. 
Debris  was  accumulated  on  the  chaimel  banks  and  lodged  in  the  vegetation. 
Debris  was  abundant  in  reaches  4  and  5,  which  is  evident  by  the  large  percentage 
differences  for  these  reaches.  The  Manning’s  n  values  for  reaches  4  and  5  in 
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I  ,1 

Table  8  I 

Percentage  Difference  Between  Manning’s  n  Values  for  November  I 
1995  and  June  1996  for  Nine  Flows  and  Five  Channel  Reaches  I 

Row  No. 

Reach  | 

J _  la 

3  3  1  1 

1 

-15 

25 

55 

60 

94  H 

2 

38 

8 

71 

93 

158  I 

3 

6 

2 

69 

57 

22  1 

4 

11 

21 

54 

66 

68  1 

5 

-1 

1 

55 

68 

109 

6 

9 

19 

54 

60 

71 

7 

5 

-6 

40 

76 

107 

8 

-18 

28 

18 

32 

44 

9 

* 

-1 

-20 

15 

26 

13 

i  Average  |  4  9  48 

60 

76 

November  were  60  and  76  percent  greater,  respectively,  than  June  due  to  the 
presence  of  debris. 


Evaluation  of  drag  coefficient 


Fischenich  (1996)  presented  Equation  2  for  estimating  Manning’s  n  for  chan¬ 
nels  in  which  the  vegetation  is  not  submersed. 


n  =  k^R 


2B 


C,Veg, 


2g 


\m 


Equation  2  is  limited  because  there  is  little  information  available  concerning  die 
drag  characteristics  of  vegetation.  Through  dimensional  analysis,  it  can  be  dem¬ 
onstrated  that  the  coefficient  of  drag  is  a  function  of  the  Reynolds  number  Re 
(Vogel  1981) 

=  =  (12) 


where 

a  =  coefficient 
I  =  characteristic  length 
V  =  kinematic  viscosity 
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The  characteristic  length  I  is  expressed  as  the  hydraulic  radius  R  in  open  channel 
flow.  However,  for  submersed  objects  Re  is  typically  based  on  the  greatest 
length  of  the  object  in  the  direction  of  flow  (Vogel  1981).  The  characteristic 
length  that  should  be  used  to  define  the  Re  is  imclear  for  an  open  channel  with 
submersed  objects  such  as  vegetation.  Because  the  local  flow  velocity  through 
vegetation  can  vary  substantially,  the  velocity  that  should  be  used  to  define  Re  is 
equally  ill  defined.  In  this  study,  the  hydraulic  radius  R  and  the  mean  velocity  V 
were  used  to  define  Re. 

The  value  of  Q  for  the  vegetation  in  the  test  chaimel  was  calculated  by  solv¬ 
ing  Equation  2  for  Q  using  the  values  of  Veg^ ,  Maiming’ s  n,  and  hydraulic  ra¬ 
dius  R  for  November  and  June  flow  tests.  The  drag  coefficient  C^was  calculated 
for  the  nine  test  discharge  rates  for  a  total  of  18  values.  The  calculated  C^was 
plotted  against  both  the  product  of  the  average  velocity  times  average  hydraulic 
radius  VR  (Figures  9  and  10)  and  the  Re  (Figures  1 1  and  12)  for  the  chaimel  con¬ 
ditions  in  November  1995  (debris  present  and  green  foliage  absent)  and  June 
1996  (debris  absent  and  green  foliage  present),  respectively.  The  plots  indicate 
that  the  Q  approaches  a  constant  value  at  high  velocities. 


Rgure  9.  as  a  function  of  VR  for  nine  test  discharges  based  on  measure¬ 
ments  conducted  in  November  1995  with  debris  present  and  leaves 
absent 
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Rgure  12.  as  a  function  of  Re  for  nine  test  discharges  based  on  measure¬ 
ments  conducted  in  June  1996  without  debris  and  with  leaves 

Curve  fits  for  the  data  were  obtained  by  regression  using  the  method  of  least 
squares.  The  relation  between  C^and  VR  under  conditions  with  debris  present 
and  leaves  absent  from  the  trees  and  shrubs  was 

Cj  =  2.1(VR)-‘’  (13) 


and  with  debris  removed  and  leaves  present  was 

Q  -  0.28(W?)-‘-*  (14) 


The  coefficients  of  determination  for  the  regression  functions  were  0.99  and 
0.98  for  Equations  13  and  14,  respectively.  A  graph  of  Equations  13  and  14  on 
log-log  scale  (Figure  13)  depicts  two  curves  having  the  same  slope;  however,  the 
presence  of  debris  in  the  charmel  resulted  in  a  drag  coefficient  C^that  was  7.5 
times  greater  than  without  debris  for  a  given  VR  product 

The  relation  between  and  Re  with  debris  present  and  leaves  absent  was 
(Figure  11) 


Q  =  9.3jcl0^(/fc)'‘-‘ 


(15) 
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Rgure  13.  C^asa  function  of  VR  plotted  on  log-log  scale  for  conditions  witfi 

debris  and  without  leaves  (Equation  13)  eind  without  debris  and  with 
leaves  (Equation  14) 

and  without  debris  and  leaves  present  (Figure  12) 


Cj  =  9.1x10^ 


(16) 


The  kinematic  viscosity  vof  the  water  was  estimated  to  be  9.1  x  lO'^  m^s  based 
on  a  water  temperature  of  13.3  “C  in  November  1995  and  1.2  x  10‘^m^s  based 
on  water  temperature  of  24.4  °C  in  June  1996.  A  graph  of  Equations  15  and  16 
on  log-log  scale  (Figure  14)  depicts  two  curves  having  the  ^proximately  the 
same  slope.  The  magnimde  of  C^for  a  specific  Re  was  10.2  times  greater  in 
November  1995  than  in  June  1996  presumably  because  of  die  presence  of  debris 
and,  to  a  lesser  extent,  the  exclusion  of  the  grass  from  the  Veg^term. 

Assumptions  of  the  methodology 

Equation  2  was  derived  from  the  de  Saint  Venant  and  Manning  equations  for 
open  chaimel  flow  based  on  several  assumptions: 

a.  Flow  depth  is  less  than  the  height  of  the  vegetation. 

b.  Vegetation  is  randomly  distributed  within  the  control  volume. 
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Rgure  1 4.  Cj  as  a  function  of  Re  plotted  on  log-log  scale  for  conditions  with 

debris  and  without  leaves  (Equation  1 5)  and  without  debris  and  with 
leaves  (Equation  16) 

c.  Resistance  due  to  surface  friction,  flow  obstructions,  and  changes  in  chan¬ 
nel  geometry  is  negligible  compared  with  vegetative  effects. 

d.  Pressure  distribution  is  hydrostatic. 

e.  Flow  velocity  over  the  channel  cross  section  is  uniform 

f.  Steady,  uniform  flow  conditions  exist. 

The  assumptions  inherent  in  Equation  2  were  employed  in  order  to  reduce  the 
number  of  required  input  parameters  and  provide  a  practical  procedure  for  ob¬ 
taining  reasonable  estimates  of  flow  resistance  and  chaimel  capacity.  Many  of 
these  simplifying  assumptions  were  violated.  Grasses,  vines,  and  shrubs  were 
sometimes  submersed  near  the  center  of  the  V-shaped  test  channel  where  the 
flow  velocity  and  depth  were  greatest  Flippin-Dudley  et  al.  (in  preparation) 
demonstrated  that  the  vegetation  in  the  test  channel  was  probably  not  randomly 
distributed  and  is  rarely  completely  random  in  nature  (Bonham  1989).  Surface 
friction  may  have  accounted  for  ^proximately  19  percent  on  average  of  the  total 
resistance  in  November  and  29  percent  in  June,  assuming  resistance  factors  are 
additive  and  a  surface  Maiming’ s  n  value  of  0.025  (Arcement  and  Schneider 
1989).  Extreme  variation  in  transverse  flow  velocity  was  observed,  as  evidenced 
by  the  high  values  of  the  energy  coefficients  presented  in  Tables  5  and  6. 

Another  assumption  of  this  study  was  that  the  coefficient  of  drag  and  density 
of  vegetation  remained  constant  throughout  the  duration  of  the  test  An  average 
vegetation  density  was  used  to  represent  the  entire  length  of  the  channel  under 
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investigation.  Both  the  vegetation  density  and  the  drag  coefficient  decreases  as 
velocity  increases  and  the  vegetative  elements  are  deflected  in  the  direction  of 
flow. 

The  combined  effect  of  the  assumptions  of  the  method  must  be  absorbed  by 
the  drag  coefficient  Laboratory  research  has  demonstrated  that  the  coefficient  of 
drag  for  flow  around  cylinders  and  similar  objects  is  approximately  1.0  (e.g., 
Petryk  1969  or  Tseng  1974).  The  drag  coefficient  associated  with  the  June 
measurements  was  approximately  1.0  for  Re  >  2  x  10^  (Figure  11),  which  indi¬ 
cates  that  the  assumptions  made  with  respect  to  the  flow  and  vegetative  condi¬ 
tions  were  reasonable  for  the  analysis  conducted  using  June  data.  The  coefficient 
of  drag  was  approximately  6.5  for  >  2  x  10®  measured  in  November.  The 
high  coefficient  of  drag  is  primarily  due  to  the  presence  of  debris  and,  to  a  lesser 
extent,  grass  that  was  not  included  in  the  parameter  used  in  the  development 
of  Equation  13  and  Equation  15. 

Evaluation  of  debris 

Flow  resistance  due  to  the  presence  of  debris  can  be  evaluated  by  extending 
Equation  2  to  include  the  effects  of  debris.  Assuming  that  the  drag  characteris¬ 
tics  of  debris  are  the  same  as  those  for  the  associated  vegetation. 


n  = 


C/Veg,*Deb,) 


11/2 


2g 


(17) 


where  Debj  is  the  frontal  area  of  the  debris  projected  onto  a  plane  perpendicular 
to  the  direction  of  flow  and  the  other  parameters  are  as  defined  previously. 

To  assess  the  value  of  Debj ,  the  drag  coefficient  Cj  for  Re>2x  10®  was 
assumed  to  equal  unity,  and  the  density  of  the  grass  in  November  was  equal  to 
that  measured  in  June.  The  Veg^  including  the  grass  was  0.322/m  ((68  grass  hits 
+  21  stem  hits)  divided  by  906  points/0.3  m).  Equation  17  was  solved  for  Deb^ 
based  on  flow  measurements  conducted  in  November.  The  average  Debris  esti¬ 
mated  to  be  0.88/m,  which  would  account  for  approximately  73  percent  of  the 
Veg^  and  Deb^  combined. 

Application  of  the  Methodology 

The  proposed  equations  for  ,  Equations  13  through  16,  can  be  applied  in 
conjunction  with  Equation  2  to  determine  the  Manning’s  n  for  chaimels  having 
vegetative  and  debris  conditions  similar  to  those  tested  in  this  study.  Arcement 
and  Schneider  (1989)  presented  data  for  several  wooded  forest  floodplains  in 
Mississippi  and  Alabama.  The  data  presented  for  Pea  Creek  (Ming,  Colson,  and 
Arcement  1979),  Coldwater  River  (Colson,  Ming,  and  Arcement  1979a),  Thomp¬ 
son  Creek  (Colson,  Arcement,  and  Ming  1979),  and  Yockanoonkany  River 
(Colson,  Ming,  and  Arcement  1979b)  were  used  to  demonstrate  use  of  the 


28 


Chapter  2  Field  Test  and  Application 


equations  to  predict  Manning’s  n  for  a  vegetated  channel.  A  sununary  of  the 
relevant  parameters  is  presented  in  Table  9. 

Arcement  and  Schneider  (1989)  presented  photographs  of  each  stream  accom¬ 
panied  by  a  description  of  the  floodplain  and  vegetation  characteristics  including 
the  Vggj,  type  of  vegetation  and  influence  of  ground  cover  and  obstmctions.  The 
flow  measurements  were  conducted  in  winter  when  no  leaves  were  present  on  the 
trees,  and  Veg^  was  determined  from  direct  measurement  of  tree  and  large  vine 
diameters.  Ground  cover  and  undergrowth,  such  as  shrubs  and  grass,  and  vegeta¬ 
tive  debris  were  not  included  in  the  Vieg^term. 

The  graph  of  the  predicted  versus  measured  Manning’s  n  values  using  Equa¬ 
tions  13  and  14  to  predict  the  drag  coefficient  C^is  shown  in  Figure  15.  The  n 
values  predicted  using  Equation  13  plot  close  to  the  line  of  perfect  agreement 
while  Equation  14  substantially  underestimated  the  n  value.  This  result  is  not 
suiprising  because  the  contribution  of  vegetative  debris  and  undergrowth  were 
not  included  in  the  Veg^  measurements  conducted  by  Arcement  and  Schneider 
(1989). 
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Rgure  1 5.  Predicted  versus  measured  Manning’s  n  for  stream  data  presented 
by  Acrement  and  Schneider  (1989) 


Chapter  2  Field  Test  and  Application 


29 


3  Conclusions  and 
Recommendations  for 
Future  Research 


Prediction  of  Flow  Resistance 
in  Vegetated  Open  Channel 
and  Floodplains 

A  vegetated  channel  located  at  tiie  ARS  Outdoor  Hydraulic  Laboratory  near 
Stillwater,  OK,  was  used  to  demonstrate  the  application  of  the  horizontal  point 
frame  and  its  role  in  a  comprehensive  approach  to  predicting  resistance  to  flow. 
The  vegetation  measurements  were  used  in  conjunction  with  flow  measurements 
in  the  test  channel  to  develop  a  relationship  for  die  drag  coefficient  of  the  vegeta¬ 
tion  Cj .  Flow  tests  were  conducted  for  two  vegetative  conditions:  (a)  in  Novem¬ 
ber  1995  when  leaves  on  the  trees  and  shrubs  were  absent  and  debris  was  pres¬ 
ent,  and  (b)  in  June  1996  when  leaves  were  present  and  the  debris  was  removed. 
An  evaluation  of  the  data  indicated  the  following: 

•  The  value  of  Cj  can  be  related  to  the  product  of  the  velocity  and  hydraulic 
radius  VR  and/or  the  Reynolds  number  Re. 

•  With  leaves  absent  and  debris  present,  Q=  2.1(VR)'‘  ' ,  Equation  13. 

•  With  leaves  present  and  debris  absent,  Cj=  0.28(VR)'*  * ,  Equation  14. 

•  With  leaves  2d>sent  and  debris  present,  Q=  9.3  x  10*(Re)'*  * ,  Equation  15, 
based  on  a  kinematic  viscosity  of  9.1  x  10"’  m^/s. 

•  With  leaves  and  debris  absent,  Cj=  9.1  x  10*(/?e)'*  * ,  Equation  16,  based 
on  a  kinematic  viscosity  of  1.2  x  lO"*  m^/s. 

•  Cj  was  approximately  unity  for  June  measurements  at  Re  >  2  x  10*  which 
is  consistent  with  laboratory  model  tests  performed  by  others. 
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•  Q  was  approximately  6.5  for  November  measurements  at  /Je  >  2  x  10®, 
greater  than  the  expected  value  of  unity  primarily  due  to  the  presence  of 
debris. 

•  Equations  13  or  15  can  be  used  in  conjunction  with  Fischenich’s  (1996) 
equation  (Equation  2)  to  determine  Maiming’s  n  when  debris  is  not  present 
and  vegetation  is  not  submersed. 

•  Equations  14  or  16  can  be  used  in  conjunction  with  Fischenich’s  (1996) 
equation  to  determine  Manning’s  n  when  debris  is  present,  leaves  on  trees 
and  shrubs  are  absent,  and  minimal  ground  cover  is  present 

•  Using  Equations  14  or  16,  the  predicted  Manning’s  n  was  within  an  aver¬ 
age  of  10  percent  of  the  measured  value  for  four  streams  in  Mississippi 
and  Alabama. 

Recommendations  for  Future  Research 


Vegetation  plays  a  vital  role  in  the  environmental  balance  of  riparian  ecosys¬ 
tems.  Standard  methods  for  characterizing  vegetation  and  predicting  flow  resis¬ 
tance  allow  interdisciplinary  teams  to  address  practical  problems  using  compara¬ 
ble  data  sets  and  analysis.  Engineers  could  collaborate  with  other  scientists  to 
determine  the  most  beneficial  plan  to  mitigate  flood  hazards  and  restore  or  main¬ 
tain  the  equilibrium  of  a  stream  systenL  Once  a  common  foundation  and  inter¬ 
pretation  of  the  design  parameters  can  be  reached,  die  findings  by  engineering 
and  scientific  disciplines  can  be  better  incorporated  into  river  management  prac¬ 
tices.  Pursuant  to  diis  end,  future  research  should  focus  on  the  following: 

•  The  use  of  Veg^  as  a  parameter  for  classification  and  optimization  of  the 
flow  resistance,  erosion  and  sedimentation,  bank  stability,  habitat  and 
aesthetic  properties  of  riparian  vegetation. 

•  The  optimum  quadrat  size  and  shape  for  the  direct  measurement  of  tree 
trunks  used  to  determine  Veg^ . 

•  The  relation  between  the  Vegj  measured  without  flow  and  the  actual  Veg^ 
with  flow  for  flexible  vegetation. 

•  Quantitative  methods  for  the  measurement  of  vegetative  debris. 

•  A  vegetative  parameter  that  characterizes  the  spatial  patterns  and  distribu¬ 
tion  of  vegetation  (i.e.,  clustering,  etc.). 
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